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Abstract 

Small perturbations of the homogeneous cooling state (HCS) for a low density granular gas 
are described by means of the linearized Boltzmann equation. The spectrum of the generator 
for this dynamics is shown to contain points corresponding to hydrodynamic excitations. The 
corresponding eigenvectors and eigenvalues are calculated to Navier-Stokes order and shown to 
agree with those obtained by the Chapman-Enskog method. The conditions for the hydrodynamic 
excitations to dominate all other excitations are discussed. 
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I. INTRODUCTION 



The use of hydro dynamic equations to describe granular fluids in rapid flow has been in 
practice for many years The justification for this fluid-like description and prediction of 
the transport coefficients appearing in these equations has been the focus of attention for 
some time as well 3, 0, Q]. In recent years, an accurate derivation of Navier-Stokes order 
hydrodynamics has been given from the Boltzmann equation for granular gases using an 
adaptation of the Chapman- Enskog method for normal gases 0, O]. The expressions for 
the transport coefficients as a function of the degree of inelasticity have been confirmed by 
both Monte Carlo and molecular dynamics simulation M. Successful application of these 
Navier-Stokes equations to a number of different states also supports their validity |7|, |8|, |9|| . 
However, the context of the hydrodynamic equations remains uncertain. What are the 
relevant space and time scales? How much inelasticity can be described in this way? 

Such questions can be addressed for gases using the Boltzmann kinetic equation to de- 
scribe the complete dynamics for properties of interest. Then it can be asked under what 
conditions do the hydrodynamic excitations emerge as the dominant dynamics . The analysis 
of this problem for normal gases is quite complete and precise Q, [ill. The objective here 
is to initiate a similar formulation of the problem for granular gases First, solutions 
to the Boltzmann equation are considered for states that deviate from spatial homogeneity 
only by small perturbations. The dynamics in this case is governed by the linear Boltzmann 
operator for spatially inhomogeneous states. The spectrum of this operator determines all 
possible excitations on all space and time scales, and for all degrees of inelasticity. The 
first problem is to identify points in this spectrum corresponding to hydrodynamics. This is 
one of the main results reported here. Both the hydrodynamic eigenvalues and eigenfunc- 
tions are calculated for long wavelength excitations corresponding to Navier-Stokes order 
hydrodynamics. Their agreement with the corresponding results from the Chapman-Enskog 
method is established. 

Next, the issue of conditions for the dominance of the hydrodynamic excitations, or 
modes, is considered. This entails showing that the hydrodynamic eigenvalues are smaller 
than all other parts of the spectrum, such that there is a long enough time scale for the latter 
to decay relative to the hydrodynamic modes. There are two new difficulties for granular 
gases. First, the hydrodynamic eigenvalues cannot be made arbitrarily small since they do 
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not all scale with the inverse wavelength of the perturbation, as for normal gases. Second, 
there is a new characteristic frequency, the cooling rate, in addition to the collision frequency 
to set the scale of the spectrum. The cooling rate can be relatively large or small depending 
on the inelasticity of the gas particle collisions. 

Analysis of the non-hydrodynamic spectrum of the Boltzmann equation for a granular 

H remains a difficult unsolved problem. Instead, we consider here a simpler kinetic model 
that retains the exact hydrodynamic spectrum and allows identification of the entire 
spectrum as well. It is found that the time scale for the hydrodynamic excitations is longer 
than that for all other excitations, for any degree of inelasticity. A brief report of these 
results has been given in reference [l^], with the details and elaboration given below. 

The plan of the paper is as follows. In the next section, a short summary of the Boltzmann 
equation is given, and the Navier-Stokes equations derived from it by the Chapman-Enskog 
procedure are recalled. More details are given in Appendix ^ Also, the hydrodynamic 
equations are linearized and the hydrodynamic modes to second order in the wavevector are 
identified. 

In Sec. mil the structure of the linear Boltzmann equation is discussed, emphasizing 
the relevance of the spectrum of the operator generating the linear dynamics to determine 
the existence and validity of hydrodynamics. The first of the above issues is addressed in 
Sec. IIV| where the kinetic theory hydrodynamic modes are identified in the long wavelength 
limit. A technical point associated with the non Hermitian character of the operator is 
discussed in Appendix |Bl The results are extended to Navier-stokes order in Sec. El and 
the obtained expressions are shown to agree with those derived in Sec. Ofrom the hydrody- 
namic equations. Details of the calculations are given in Appendix O The possibility of a 
description in terms of the hydrodynamic modes is studied in Sec. IVIl by means of a model 
kinetic equation. It is established that there are length and time scales on which only the 
excitations associated to the hydrodynamic modes persist. The mathematical details are 
given in Appendix O Finally, the last section contains a short summary of the main points 
and conclusions in the paper. 
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II. BOLTZMANN EQUATION AND HYDRODYNAMIC MODES 



In all of the following, the simplest model of a granular fluid is considered: a low density 
gas of smooth and inelastic hard spheres {d = 3) or disks {d = 2) of mass m and diameter a. 
The binary collision rule given below is characterized by a constant coefficient of restitution 
a, defined in the interval < a < 1, and measuring the loss of energy in each collision. At 
sufficientlvlow density, the distribution function /(r, v, t) is determined from the Boltzmann 
equation Q, 14 , 3| 

{dt + v-V)f = J[fJl (1) 
where J is the inelastic Boltzmann collision operator defined by 

J[X, Y] = a"-' j dv, j da Q{a 9)^ 9 [a-^X{r, v', t, )Y{r, v[, t) - X(r, t)F(r, v^, t)] , 

(2) 

for arbitrary functions X{v) and Y{v). Here, ct is a unit vector along the line joining the 
centers of the colliding pair, O is the Heaviside step function, and g = v — Vi. The primes on 
the velocities denote the initial values {v[, V2} that lead to {vi, V2} following a "restituting" 
binary collision. 



v' = v -]-{! + a ^)(ct ■ g)a, v[ = Vi + hi +a ■ g)a. 



(3) 



The usual Boltmann collision operator is recovered from Eq. ^ in the elastic limit a = 1. 

The macroscopic variables of interest are the hydrodynamic fields, i.e. the density n{r, t), 
the fiow velocity u{r,t), and the (granular) temperature T{r,t). They are defined as mo- 
ments of the solution to the Boltzmann equation. 



n(r, t) 
n{r, t)u(r, t) 



dv 



1 

V 



f{r,v,t). 



(4) 



yin{r,t)T{r,t) J 

An exact set of equations for these variables is obtained from the following properties of the 
collision operator 



dv 



1 

V 



y |m [v — u) 



( 



J[f, f] 








\ 



\-lnTC) 



(5) 



The first two equations follow from conservation of mass and momentum in the particle 
collisions. The last equation reflects the loss of energy in collisions due to the inelasticity. 
This appears through the "cooling rate" C(^j^) ^ 0, defined by this equation. 
The macroscopic balance equations resulting from the above properties are: 

dtn + V ■ (nu) = 0, (6) 

{dt + u-V)u+ (mn)-^ [V{nT) + V ■ □] = 0, (7) 
2 

(dt + u-V + C)T + — (riTV • u + □ : Vu + V ■ g) = 0. (8) 

nd 

The functionals giving the dissipative part of the pressure tensor, fl, and the heat flux, q, 
are also moments of the solution to the Boltzmann equation, 



n(r,t) = Jdvm(vV-^\^ fir,v,t), 



(9) 



m ^ ^2 d + 2^ 



q{r, t) = J dv \^-V' - -^Tj V f{r, v, t), (10) 

where I is the unit tensor of dimension d and V = v — u the so-called peculiar velocity. 
Equations (jE))-® are the basis for a hydrodynamic description, once n(r,t), q{r,t), and 
C{r,t) are specified. These can be obtained from their definitions using a solution of the 
Boltzmann equation generated by the Chapman- Enskog method. This method assumes 
the existence of a solution whose space and time dependence is given entirely through the 
hydrodynamic fields and their gradients (i.e., a "normal" solution). As a result, n(r,t) 
and q{r,t) are given in terms of these variables, and Eqs. (jEI)-® become a closed set of 
hydrodynamic equations. The primary results of this method are recalled in Appendix 1X1 
To leading order in the spatial gradients, the dissipative fluxes are found to be given by 

m 

njj(r, t) = -T] {ViUj + VjUi - ^(^ij-V ■ 11^ , (11) 
q(r, t) = -kVT - /iVn. (12) 

For the cooling rate the result is: 
/2T\ 

C(r, t) = na'^"^ ( — j Co + CiV^T + C2V^n + bilinear in Vn, VT, Vu terms, (13) 



5 



with Q being a dimensionless positive constant proportional to (1 — a^). The exphcit forms 
for the shear viscosity 1], the thermal conductivity k, and the new transport coefficients /i, 
(i, and (2, peculiar to granular gases, are given in Appendix^ The nonlinear contributions 
to indicated in (jl3|) play no role in the following linear analysis and will be not discussed 
further in this paper. 

Equations (0)-® together with the "constitutive relations" (fTT |) - (fT^ are the Navier- 
Stokes hydrodynamic equations for a granular gas. A special solution of these equations is 
that for a homogeneous state, characterized by a constant density rihcs, a vanishing velocity 
flow Uhcs = 0, and a uniform temperature Thcsif) determined from 

dtn^s it) = -n,^sa'-' (1^) CoT',i; it) . (14) 

This is referred to as the homogeneous cooling state (HCS), since the only macroscopic 
dynamics is the monotonic decrease in the temperature field. It is easily seen that the HCS 
as defined above is also a solution of the exact balance equations, taking into account that 
(hcs oc nhcsO'^'^iThcs/nT')^^'^ on dimensional grounds. In the following, the solution to the 
Boltzmann equation will be considered for small perturbations of the HCS, and it will be 
useful to have the corresponding results from Navier-Stokes hydrodynamics. The linearized 
Navier-Stokes equations have time dependent coefficients, since the distribution function of 
the reference HCS depends on time. This complication can be overcome by the introduction 
of suitable dimensionless space and time scales, as well as scaled hydrodynamic fields. Thus 
we define r* and s by 

r* = j, ds=^-^dt. (15) 

Here, Vhcs it) = [2Thcsit)/'m]^^'^ is the thermal velocity in the HCS, and £ = i'rihcs'^'^^^)^^ is 
proportional to the mean free path of the particles. The dimensionless fields 6yjir*, s) are 
chosen as 

/ * N nir,t)-nhcs * ^ T(r, t) - T/,cs(t) , , . uir,t) 
p{r , s) = , eir , s) = — — , u:ir ,s) = --. (16) 

nhcs Thcs{t) Vhcs{t) 

Besides, since the equations are linear and the reference state is homogeneous, it is sufficient 
to consider a single Fourier mode, 

6y,ir\s) = e^''-^''6y,ik,s). (17) 

The velocity components are chosen as a longitudinal component relative to k, uj\\ (fc, s) = 
k uj (fc, s), and d—1 transverse components u]±^iik, s) = e*-a3(A;, s), where {k = k/k, e^, i = 



1, . . . ,d—l} are a set of d pairwise orthogonal unit vectors . The dimensionless hnear Navier- 
Stokes equations then become a system of ordinary differential equations with constant 
coefficients that can be expressed in the compact form ^ 

dJy{k,s) + K{k)-Sy{k,s) = 0, (18) 

where we use a + 2 dimensional space representation with 6y{k, s) being the vector 

fp{k,s) \ 

e{k, s) 

uj\\{k, s) 
\u}±{k,s) j 

and C(3_L(fc, s) denoting the vector formed by the d — \ components tD_L^j(fc, s). The matrix 
K (/c) is block diagonal, expressing the decoupling of transverse and longitudinal modes, 



by{k,s) 



(19) 




Ki(fc) 



v 



K(A:) 





Co + - C2) f + - cr) 



(20) 



2 ^ 



ik 
^k 

a 

2(q!-1) ^*l2 
d 



^0 



Uk) = -[^-^-r^*k']\. 



(21) 



(22) 



In the above expressions, K2 and I are matrices of dimension d — 1. The dimensionless 
transport coefficients used in the above equations are defined by 

V * 2fi; 2/i 



V 



c 



K 



diUhcsVh c 



/i 



diThcsVhc 



(23) 



cr 



ThcsC 



1 

~5 S2 



nhcsC2 



^"Vhcs ^"Vhcs 

The formal solution to the initial value problem ()18p is 



dy{k,s) 



-K(fc)s 



6y{k,0). 



(24) 



(25) 



The eigenvalues and eigenvectors of the generator K for this dynamics define the d + 2 
Navier-Stokes order hydrodynamic modes. They are given by the solutions of the equation 



K(A;) ■ <p,{k) = Xj{k)ip,{k), j = l,...,d + 2. 



(26) 



A simple calculation provides the expressions for the eigenvalues to order k"^. They are given 

— Co fd + l , ,A 



Ai(A:) = -, A2(A;) = ^ 



A||(fc) 



^0 



+ 



1 I 2(ci-l)r7- 



Co 



(27) 



the eigenvalue Ax(A;) being (rf— l)-fold degenerate. The corresponding eigenvectors to leading 
order in k are 

/ 1 \ 

-2 





V / 



/o\ 

1 



/o\ 






(^11 (fc) 





1 



(28) 



Here, 



= 0, and i = 1 = 1, for = 2. 





and 2 



for (i = 3. 



The first of these modes is excited by the condition 9 {k, 0) = —2p{k, 0) at zero flow 
velocity. The second is produced by a temperature perturbation at constant density and 
also zero velocity, while the third one corresponds to a longitudinal velocity perturbation at 
constant temperature and density. There is a (rf — l)-fold degeneracy for the shear modes 
of eigenvalue A_l. These diffusive modes are excited by a perturbation of the velocity field 
in the transverse plane orthogonal to k. 

It should be noted that while the above analysis is restricted to the Navier-Stokes equa- 
tions, derived from the Boltzmann equation by the Chapman- Enskog method, the eigenval- 
ues and eigenvectors to order k follow more generally from the exact macroscopic balance 
equations and do not depend on the approximate constitutive equations (fTTjl and (fT^ . The 
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hydrodynamic modes sought by kinetic theory in the subsequent sections can therefore be 
defined as those excitations due to small perturbations which agree with the above in the 
long wavelength limit. Analyticity then allows extension of that identification to shorter 
wavelengths. A consistency check of the Chapman-Enskog method is agreement with the 
above results at order k"^. However, the concept of hydrodynamic modes in this context 
does not require the validity of the Chapman-Enskog method nor the limitation to the 
Navier-Stokes approximation. 



III. LINEAR BOLTZMANN EQUATION 

A more complete and accurate description of the response to small spatial perturbations of 
the density, temperature, and fiow velocity is obtained directly from the Boltzmann equation. 
Consider first an isolated system. As already indicated in the previous section, the exact 
balance equations have a solution describing the HCS, with a monotonically decreasing 
temperature obeying Eq. (fT^ . The solution to the Boltzmann equation corresponding to 
this macroscopic state is characterized by the scaling form jj] 

fhcs{v,t)=nj,,,v^^^^{t)(j){V*), (29) 

with 

V* = = (30 

Vhcs{t) Vhcs[t) 

Substitution into the Boltzmann equation gives 

'^^^Y ' C^fhcs) = J[fhcs, fhcs]- (31) 

For later convenience, a constant velocity Uhcs for the system as a whole has been included, 
although this can always be removed by means of a Gallilean transformation. Then, in the 
following it will be considered that u^cs = 0, unless it be explicitly established otherwise. An 
exact and explicit solution of Eq. (jHT|) is not known yet, but the behavior of (p at large and 
small velocities has been determined uM and the results obtained by the direct simulation 

n 

Monte Carlo method strongly supports the existence of such a scaling form jl6| . 

The HCS distribution function is a "universal" homogeneous solution in the same sense as 
the Maxwellian for elastic collisions. An arbitrary homogeneous state is expected to approach 
the HCS after a few collisions. Therefore, in discussing response of any homogeneous state 



to small spatial perturbations, it is sufficient to consider the HCS as the reference state. All 
the other cases will simply induce additional short time transients. 
Consider then small spatial perturbations of the HCS 

f{r,v,t) = fhcs{v,t)[l + A{r,v,t)], \A{r,v,t)\ 1. (32) 

To linear order in A(r, v,t), the Boltzmann equation becomes 

[dt + vV + L{t)] iUcsA) = 0, (33) 

where L{t) is the linearized Boltzmann collision operator given by 

L{t)X{v) = -J[fhcs,X] - J[X,Ucs], (34) 

for arbitrary X{v). Just as for the analysis of the Navier-Stokes equations in the previous 
section, the above linear kinetic equation takes a simpler form when expressed in terms of 
the dimensionless variables defined in Eq. (fT^ . and considering a single Fourier mode 

A(r, V, t) = e^*"^* A(fc, v*, s), (35) 

where v* = v/vhcsif). Then, Eq. PHj) becomes 

{ds + ik-v* + C*)'K{k,v*,s) = Q. (36) 

The dimensionless operator C* is now time independent and it is given by 

C*X = ^(j)-^^ ■ (v*(l)X) + L*X, (37) 
2 ov* 

where L* is the dimensionless linear Boltzmann collision operator 

L*X = -0-1 ( J*[0, (PX] + J*[(j)X, (f)]) , (38) 



r[x,Y] ^ j dvij dae{a-g*)a-g* [a-2x(0>^«) -^K)>^K)] • (39) 

Here, v*' and v^' are related with v* and by Eqs. (jH)). The operator C* differs from 
the linearized Boltzmann collision operator L* by terms representing the cooling effects of 
the inelastic collisions. The latter arise because the derivative with respect to s is taken at 
constant v* rather than v. 
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Solutions to the dimensionless, linear kinetic equation ()36p for A are sought in a Hilbert 
space whose scalar product is defined by 

{X,Y) = J dv*(P{v*)X^v*)Y{v*), (40) 

with the dagger denoting complex conjugation. The formal solution to the kinetic equation 
is 

A (fc, V*, s) = — Idz e-''n{z)A* (k, V*, 0) , (41) 
2111 J 

n[z) = {z-ik-v* - C*y^ , (42) 

where the contour encloses the entire spectrum of ik ■ v* + C* , both point and residual, 
counterclockwise. It is important to realize that all the linear excitations of the granular 
gas are determined from this spectrum. This formulation of the problem for small spatial 
perturbations provides a precise context for addressing many questions regarding hydrody- 
namics for a granular gas. The existence of the hydrodynamic modes and their role relative 
to other dynamical processes are determined by the characterization of the above spectrum. 
To see how this occurs, it is useful to recall briefiy the status of the corresponding problem 



for gases with elastic collisions 



IG 



IjJ. In that case, it has been proved that the hydro- 
dynamic modes exist as d + 2 poles located at the origin in the long wavelength limit and 
corresponding to the local conserved quantities. Furthermore, the remainder of the spec- 
trum is bounded away from these poles, and the spectrum is analytic in k about = 0, so 
this isolation of the hydrodynamic modes is preserved at finite wavelengths. 



IV. EXISTENCE OF HYDRODYNAMIC MODES 



The spectrum of ik ■ v* + C* is expected to be quite complex, based on the special case 
of elastic collisions, with points spectra, continua, and limit points. The hydrodynamic 
excitations, whenever they exist, are part of the point spectrum so in order to investigate 
them it suffices to consider the eigenvalue problem 

{ik ■ V* + C*) = Xi (k) (43) 

The search for hydrodynamic excitations can be carried out by assuming they are analytic 
in k and looking first for the fc = solutions of Eq. ()43|) . The practical issue of constructing 
these modes at finite k is addressed in the next section. 
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The central idea for constructing the hydrodynamic eigenvalues and eigenvectors at k = 
is the note that the HCS is parameterized by the hydrodynamic fields Uhcs, Thcs, and Uhcs, 
which is now considered different from zero, as discussed at the beginning of the previous 
section. Therefore, differentiating the Boltzmann equation for the distribution function 
of the HCS, Eq. (jSH), with respect to these fields gives exact properties of the linearized 
Boltzmann collision operator. For example, 

d f Chcs d 



duhc 



2 dV 



i^fhcs) — J[fhcs, fl 



hcs\ 







gives directly 



■ Jhcs) + 



V 



hcs 



+ L 



df> 



hcs 



0. 



(44) 



(45) 



2dnhcsdV 2 dV \ dnh^J duhcs 

The derivatives of (hcs and fhcs are easily calculated from the properties (hcs oc Uhcs and 
fhcs oc Uhcs- In terms of the dimensionless variables, and setting Uhcs = 0, Eq. becomes 



[V 



2 ' dv* 

Next, calculate the derivative of Eq. with respect to T^^s, 



IdChcs d X , Chcs d 

■ [V Jhcs) + 



2dThcs dV ' ' ' 2 dV 
Since fhcs has the scaling form (|^. it is 

d fhcs fhcs 



V 



df^ 



hcs 



dT, 



hcs 



+ L 



dfi 



hcs 



dT, 



0. 



hcs 



d+V ■ 



d\nf 



hcs 



dThcs 2Thcs V dV 

il/2 



(46) 



(47) 



(48) 



and, taking into account that C,hcs oc ' becomes 



Chcs d 

2 dV 



(Vf 



hcs) 



Chcs d 



2 dV 



- L 



fl 



hcs 



Vfhcs [d + V 
d In fhcs 



d\nf 



hcs 



d+V- 



dV 



dV 
= 0. 



(49) 



Setting Uhc 



and transforming to dimensionless variables, the above equation yields 
din 



C* d + V* 



Co f, * din, 
^{d + v* 



dv* J 2 \ dv* 
Finally, differentiating Eq. (jHlll with respect to Uhcs, afterwards setting Uhc 
ducing dimensionless variables leads in a similar way to 



din 4 
dv* 



22. 

2 



dlncf. 
dv* 



(50) 
0, and intro- 

(51) 
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The d + 2 equations (j^ . (jSUjl . and ^T\i provide exact properties of £*. In fact, Eqs. (HI)) 
and (j^Tj) are of the form of the eigenvalue problem to be solved at = 0, with eigenvalues 
given by Co/2 and — Co/2, respectively. It is straightforward to construct linear combinations 
of (j46|) and (j5U|) to obtain an additional eigenvalue and eigenvector. The results can then 



be expressed as 



12| 



£*^i(0) = A, (0)^,(0), z = l,...,d + 2, (52) 

with 

{A,(0)} = |o,|,-|}, (53) 

{*.(0)} = |d + 1 + . - . . --^1 . (54) 

The eigenvalue — Co/2 is cZ-fold degenerated. For convenience for the perturbation calculation 
to be carried out in the following section, the subspace associated to it has been rearranged. 
The velocity v* has been decomposed into its component in the direction of fc, f jj = fc • f*, 
and the remaining d — 1 ones, forming with it a set of d pairwise orthogonal components, 
i.e. v]_ is defined by the d — 1 components v* ■ e^*). When appropriate, the set of the d 
eigenf unctions associated with the eigenvalue — Co/2 will be denoted by ^'3(0). 

Clearly, the above are the long wavelength limit of the hydrodynamic modes defined by 
Eq. ()27|) . This is a the first primary result of the analysis developed here, demonstration 
of the existence of hydrodynamic excitations in the spectrum of the linearized Boltzmann 
equation. The results are exact and apply for arbitrary degree of dissipation. 

These long wavelength eigenfunctions are determined from the HCS distribution which 
depends only on the magnitude of v*, so the terms on the right side of Eq. (|54j) are all 
determined from dln(f){v*) /dv*. In the elastic limit, (t){v*) becomes Gaussian, and the 
hydrodynamic eigenfunctions become linear combinations of 1, f *, and f i.e. of the sum- 
mational invariants for the conservation laws of mass, momentum, and energy, as expected. 
For inelastic collisions, the eigenfunctions are quite different, particularly at large v* due to 
an overpopulation in the HCS distribution relative to the Gaussian. This is illustrated in 
Figs. Hand 121 where dhicj) {y) /dv* has been obtained from Direct Simulation Monte Carlo 
solution of the Boltzmann equation for the HCS jl?! . 

The set of eigenfunctions \E'i(0) span a. d + 2 dimensional subspace of the Hilbert space, 
but they are not orthogonal. Consistently, it is easily verified that the operator C* is 
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12 3 4 
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FIG. 1: Plot of (ln0)' = dln(j)/dv* as a function of v* for d = 3 and a = 0.95. The circles are 
the numerical derivative of the DSMC results and the solid line is the Gaussian. Quantities are 
measured in the dimensionless units defined in the text. 

not Hermitian. Therefore, it is useful to introduce a biorthogonal set of functions = 
1, . . . ,d + 2, with the requirement 

i<^u^M) = k- (55) 

To identify the appropriate set of functions, first note that 

(l,rX) = 0, iv\C*X) = -^{v*,X), (56) 

as a consequence of the number of particles and momentum conservation in the moment 
equations (0). Therefore, 1 and v* are eigenfunctions of the adjoint of £*, with eigen- 
values and -Co/2, respectively. This gives the set 

{$,} = |l, $2, (57) 

The final choice of $2 does not appear to be unique. This is discussed further in Appendix 
lEl For the purposes of the next section, it suffices to make the choice 

m = !^h^ + l,k.v\vlY (58) 

The function $2 is not an eigenfunction of the adjoint operator C*^ , but it is easily verified 
that the biorthogonality conditions are satisfied. 
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FIG. 2: The same as Fig. [l] but for a = 0.6. 
V. NAVIER-STOKES ORDER MODES 

In the previous section, the hydrodynamic modes were identified in the long wavelength 
limit. Assuming analyticity, their existence at finite wavevectors can be inferred. Fur- 
thermore, their explicit construction is possible by perturbation theory. This construction 
provides a critical test of the internal consistency of other quite different approaches (e.g., 
the Chapman- Enskog method discussed above and in Appendix El- In particular, the de- 
tailed form of the eigenvalues and the dependence of the associated transport coefficients on 
the restitution coefficient should be exactly the same. This is demonstrated to Navier-Stokes 
order in this section. 

Return now to the eigenvalue problem PHj) and consider the case of k « 1. Look for 
solutions with the expansion (a more complete characterization of the conditions for this 
perturbation expansion is given below), 

^,(k) = ^,(0) + k'^f^ + P^f + . . . , (59) 

Ai(k) = A,(0) + k\f~^ + fc^Af + . . . . (60) 

The reference eigenfunctions \E'i(0) and eigenvalues Aj(0) are taken to be the long wavelength 
hydrodynamic results of Eqs. and (|^. As already indicated, there is a ci-fold degeneracy 
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for the eigenvalue A(0) = — Co/2. However, the d- dimensional subspace spanned by its 
eigenvectors was naturally partitioned by symmetry into the longitudinal and transverse 
components. The eigenvectors for the transverse modes decouple from the remaining three 
modes even at finite wavevector for the same symmetry reasons, so there are no complications 
of degenerate perturbation theory. In the longitudinal subspace, all the eigenvalues are 
distinct, except in the elastic limit a = 1, where there is a three- fold degeneracy in this 
subspace. Thus the two cases of unperturbed reference states with a = 1 and a < 1 must 
be distinguished. In the former case, the eigenvalues behave as 

\i{k, a) = U(0) + A;Af ^ + k^Xf + . . . 1 + terms of order (1 - a), (61) 

L J a=l 

and are regular in (1 — a). This occurs when the degree of dissipation because of inelasticity 
is small relative to the effects of the spatial variation. The corresponding eigenvalues are 
then similar to those of a normal gas, with d — 1 shear diffusion modes, two sound modes, 
and a heat diffusion mode. Here attention is restricted to the more interesting and relevant 
second case of fixed a < 1 with small spatial perturbations. It will be seen that the modes 
are now qualitatively different, consistently with the results reported in Sec. |Hl since the 
degeneracy of the elastic limit is lifted at the outset by the finite dissipation. 

To set up the perturbation expansion, projection operators Vi for the biorthogonal set 
{<l>i, \&i(0)} are defined by 

V,X = ^,{0) X), (62) 

for an arbitrary element X in the Hilbert space. The eigenvalue problem ()43p can be 
rearranged as 

[C* - A, (0)] (fc) = [A, (k) - A, (0) - tk ■ V*] (fc) . (63) 

Then, operating on both sides of this equation with Qi = 1 — Vi, and using the property 
[C* - Xi (0)] P^ = 0, this becomes 

[C* - A, (0)] Q,^, (fc) = [A, (k) - A, (0) - ■ v*] (fc) . (64) 

By construction, the right side of the above equation is orthogonal to the null space for 
the adjoint of [C* — A, (0)] , and the Fredholm alternative assures solutions to this equation 
jisl ]. The eigenvalue problem for \E'j (fc) is determined only up to an overall scale factor, 
amounting to the choice of normalization. It is convenient to choose 

($„vl/,(k)) = l, (65) 
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implying j = for n > 1. With this, Eq. (jMj) gives two sets of equations for the 

eigenvalues and eigenvectors, 

A, (fc) = A, (0) + ik ■ v*<il, (k)^ + (^<l>„ [C* - A, (0)] Q,^, (k)^ , (66) 



Q,^, (fc) = {Q, [C* - A, (0)] Qi}-' Qi [A, (A;) - A, (0) - tk ■ v*] (k) . (67) 

To zeroth order in k, these equations give the hydrodynamic modes of the last section in 
the long wavelength limit, consistently. To first order in k, the eigenvectors are 

vl/f)(fc) = QM'\k) = {Q,; [C* - A, (0)]}-^ Qi (a?) - zk ■ V*) M/, (0) 

= - {Q. [C* - A, (0)]}-' Q,tk ■ v*^, (0) . (68) 

The first equality in the above transformations is a consequence of the normalization con- 
dition ()65p . For the first order eigenvalues it is found: 

Af ) = (^<l>,, ik ■ v*^, (0) ^ + C*qM'^) = 0. (69) 

Now the last equality follows from the fact that each $j and \E'j (0) have the same definite 
parity under the change v* —v* and the distribution function of the HCS, 0(f *), defining 
the scalar product is invariant under this change. Thus the first term on the right side 
vanishes. The second term also vanishes for similar reasons, since $j and \E'-^^ have opposite 
parity and C* is invariant under the change in sign of the velocity. 
To second order in k the eigenvectors and eigenvalues are given by 

Q.*r = - {Q. [C* - A, (0)]}-^ Q,zk ■ \ (70) 



Af) = (^^,,tk ■ + ($,,rQ,vi/f)) . (71) 

The above expressions for the eigenvalues are evaluated in Appendix O The results have the 
same forms as given in Eqs. (jSZj)- Furthermore, the reduced transport coefficients rj*, k*, and 
Ci are determined from the same integral equations as following from the Chapman- Enskog 
solution summarized in Appendix El This confirms that the hydrodynamic modes deter- 
mined from the spectrum of the linearized Boltzmann equation are consistently determined 
to Navier-Stokes order by both methods. 
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VI. "AGEING" TO HYDRODYNAMICS 



The existence of hydro dynamic excitations only assures that there is a hydrodynamic 
contribution to the dynamics of small perturbations of the HCS. To establish a description 
in terms of these hydrodynamic excitations alone, it is necessary to characterize the rest of 
the excitations in the spectrum. For gases with elastic collisions, it has been shown that for 
sufficiently small k the hydrodynamic excitations are smaller in magnitude than all other 
excitations, and bounded away from them jiol . Consequently, there is a time scale 

beyond which only the hydrodynamic excitations persist. It is on this space and time scales 
that hydrodynamics in the usual sense applies. Typically, the conditions are wavelengths 
larger than the mean free path and times later than the mean free time. This leaves a large 
domain of macroscopic space and time scales for hydrodynamics. 

The extension of this concept of "ageing to hydrodynamics" for granular gases is expected, 
but its verification is not so straightforward. The mathematical analysis for elastic collisions 
does not transfer to the granular gas due to the significant differences in the linear collision 
operator. There are qualitative differences in the hydrodynamic modes. For example, the 
fact that energy is not conserved means that the hydrodynamic excitations cannot be made 
arbitrarily small simply by making k small. The fastest decaying hydrodynamic modes 
is that with eigenvalue Co/2 at long wavelengths. This has its maximum value at large 
dissipation, and the question arises as to whether the time scale for this mode can become 
comparable to or exceed those of the non-hydro dynamic modes at strong dissipation. 



A. Model Kinetic Equation 

Current analysis of the spectrum of the linearized Boltzmann operator for granular gases 
appears to be limited to the hydrodynamic excitations discussed here, with no characteriza- 
tion of the rest of the spectrum as yet. Consequently, in the remainder of this presentation 
these questions are addressed in the context of a model kinetic equation . This model |3| is 
an extension of the familiar Bhatnager, Gross, Krook (BGK) single relaxation time model 
' |. The Boltzmann equation can be formally written in the form 

d 



for normal gases 



^^+v-V)/ = -z/(/-^7). (72) 
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There are two significant differences of tlie model considered fiere witli respect to tlie original 
Boltzmann equation. First, the collision frequency is replaced by a velocity independent 
function of the local density and temperature, v = v{n^T). Second, the gain term of the 
Boltzmann equation is replaced by vg, where g is taken to be a Gaussian function of the 
velocity. 



g{r,v,t) = n 



h{T) 



d/ 



2 



e-^(^)^^ (73) 



vr 

Here V is the peculiar velocity defined bellow Eq. (fTUI) . The parameter h is chosen to be 
function of the local density and temperature so as to enforce the moment conditions (0) 
above. This leads to the identification 

In this way, it is assured that the exact macroscopic balance equations (jm)-(|lUj) are pre- 
served by the model. Then the Chapman-Enskog method leads to the same Navier-Stokes 
hydrodynamic equations as for the Boltzmann equation, with only the transport coefficients 
being different. In the following, it will be considered that the expression for the collision 
frequency u is chosen as scaling with nT^^"^, in order to mimic the hard sphere behavior. 
Dimensional analysis then implies the same scaling for the cooling rate (. Moreover, note 
that consistency of the model kinetic equation requires that ( < for all values of a. A pos- 
sible choice for the cooling rate is to be the same as obtained from the Boltzmann equation 
by using a local Maxwellian for the distribution function. In the same spirit, the collision 
frequency can be fixed by fitting one of the transport coefficients of the model to that ob- 
tained from the Boltzamnn equation by the Chapman-Enskog procedure in the first Sonine 
approximation. If the shear viscosity r] is used, the above choices lead to 

, , (3 - 3a + 2rf)(l + a) 

= ^ -Mr,t), (75) 

ar,t) = ^ '-^ -Mr,t), (76) 



where i>Q{r,t) is an average local collision frequency, 

i2 + d)T{d/2) \m 



Mr,t) = ————- - . (77) 



The above expressions yield 



ar,t) ^ (2 + rf)(l-a) ^ 2 + d ^ 

iy{r,t) 3-3a + 2d - 3 + 2d ' ^ ^ 
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in agreement with the model consistency requirement. 

This kinetic model reduces to the BGK model for normal gases at a = 1 Ql- Otherwise 
it reproduces all of the qualitative features of the granular Boltzmann equation, including 
a nontrivial HCS and the same hydrodynamic excitations discussed in the sections above. 
In some respects, the model kinetic equation is more complex than the Boltzmann equation 
since the collision operator is a nonlinear functional of / through the dependence of on T 
and u. However, the linearized model kinetic equation for small perturbations of the HCS 
is considerably simpler than that for the Boltzmann equation, as shown below. 



B. Model HCS and Linear Model Kinetic Equation 



The HCS equation (p?T|l for this model becomes in dimensionless form 



Co d 



{v*<j>) + u*<p = [tt (1 - CoV^^o)] 



*\-]-d/2 



exp 



,*2 



(79) 



where Co = Khcs/vhcsif), as in the preceding sections and, consistently, Vq = d-Vhcs/vhcsif) 
The solution to the above equation is 



(^*)=^o[^(l-aK)]"'^' / ^^exp 

^0 



exp 



CoV^o 



(80) 



It is easily verified that this distribution exhibits algebraic decay for large velocities, 

/ „ ^^^^^ 



P 



r 



(81) 



2 \P-2J V 2 

with p = 2z/q/Cq. Therefore, moments of degree p or greater do not exist. 

The linearized kinetic model equation for small perturbations of the HCS is obtained in 
Appendix^ with the result, in the dimensionless variables of the previous sections. 



{ds + tk-v*+C*JA{k,v*,s) = 0. 
The linear collision operator in this case is 

i ^ ' 

The projection operators Vi are the same as defined in Eq. ()(i2|) and 

V = ^Vi, Q = l-V. 



(82) 



(83) 



(84) 
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In the above expressions, the summations are over the hydrodynamic modes. The first term 
on the right hand side of Eq. (jH^ is the projection onto the hydrodynamic eigenf unctions, 
while the second one is orthogonal to this subspace. Consequently, the spectrum of £^ for 
the model kinetic equation has the same k = hydrodynamic eigenfunctions and eigenvalues 
as the Boltzmann equation, i.e. it is 

Cl,^iiO) = K (0)^,(0), t=l,...,d+2. (85) 

Furthermore, the structure of £^ decomposes into operators defined in the hydrodynamic 
subspace and its orthogonal complement. This allows more detailed analysis of the non- 
hydrodynamic spectrum. The associated eigenfunctions lie in the orthogonal complement 
Q. Consider the general form for = Q\E^ 

(*g, [c - J2 mn'^o) = (*Q, [•'i + f ■ fv) *q) 

The second term in the scalar product above can be simplified as 

|(...^*,.*,).|(i....^*^) 

= |(4'«,*Q)-|(r'9;,9'a, (87) 

where 

9oi^l=nj;yit)gUv,t)>0. (88) 
Then Eq. (jHBj) gives the desired inequality 

(vi/Q,[£;,-E.A.(o)P.]vi/Q) _ u* 

(^q,*q) 2 
Therefore, the non-hydrodynamic spectrum of £^ consists of points (or continuum) with 
real parts larger than i'q/2. The fastest decaying hydrodynamic eigenvalue is at Q/2. The 
hydrodynamic excitations are isolated from the rest of the spectrum for Q < Uq, as it is 
always the case. Assuming analyticity in k, the hydrodynamic modes at finite wavelength 
also will be isolated from the rest of the spectrum for sufficiently small k. 



1 + 



(^c 



(^Q,^C 



> 



(89) 
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In summary, the kinetic model considered here illustrates the expected behavior for the 
Boltzmann equation. The linearized kinetic equation for small perturbations of the HCS 
characterizes the complete complex response. Among the excitations there are d + 2 hy- 
drodynamic modes. At long wavelengths, these modes are isolated from and have smaller 
eigenvalues than the rest of the spectrum. Hence, there is a sufficiently long time scale on 
which only the hydrodynamic excitations persist. 

This kinetic model can be used to explore the role of hydrodynamics in great detail 



as the linear equation can be solved exactly j2l|. For example, it has been shown that 
the hydrodynamic modes extend to very short wavelengths far beyond the validity of the 
Navier-Stokes approximation. 

VII. DISCUSSION 

The objective here has been to explore the role of hydrodynamics for a granular gas by 
a direct analysis of the spectrum of the linear inelastic Boltzmann equation. This analysis 
has been shown to lead to results equivalent to those obtained previously based on the 
Chapman-Enskog method to solve the Boltzmann equation. However, the current method 
is more straightforward and less susceptible to subjective questions about the applicability of 
the Chapman-Enskog method to granular gases. In addition, the formulation of the problem 
in terms of the spectrum of the linearized collision operator is the proper setting to explore 
the context for dominance of a hydrodynamic description. 

The eigenvalue problem posed in Sec. IIVI together with the macroscopic balance equa- 
tions leads to a precise definition of the hydrodynamic modes. Here, it has been established 
that they exist for sufficiently long wavelengths. It remains to explore the form and extent to 
which they are meaningful at shorter wavelengths, and other methods to solve the eigenvalue 
problem are available to complement the simple perturbation theory applied here. This has 



been done in ref. 2l| for the kinetic model introduced in Sec. |3 demonstrating the exten- 
sion of the hydrodynamic modes to wavelengths an order of magnitude shorter than those 
required for the Navier-Stokes approximation . 

A second relevant result of the analysis here is the identification of the hydrodynamic 
eigenfunctions. This allows calculation of the hydrodynamic component for properties of 
interest, and provides the means to explore dynamical mechanisms beyond the Boltzmann 
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description based on hydrodynamics (e.g., mode coupling phenomena for fluctuations). For 
normal gases, these eigenfunctions are the summational invariants (1, f ^, v) in the long 
wavelength limit. For a granular gas, they are replaced by derivatives of the logarithm of the 
HCS distribution function, which behave quite differently for large velocities, as illustrated 
in Figs, m and |21 

The hydrodynamic excitations are usually interpreted as a set of c? + 2 modes associated 
with the density, temperature, and flow field of the macroscopic balance equations. Here, 
they are simply d + 2 points associated with the spectrum of the linearized inelastic Boltz- 
mann equation for response to small perturbations. It is often claimed that the temperature 
should not be included in the set of hydrodynamic fields, as the energy is not conserved. 
However, the analysis here shows that in the proper reduced variables all but one eigenvalue 
is non-zero in the long wavelength limit. Thus, the exclusion of one field from the descrip- 
tion will not recover a smaller set of eigenvalues clustered around zero at long wavelengths. 
Instead, the d + 2 modes have the following qualitative behavior. They are all clustered 
near zero for weak dissipation, but become increasingly separated at stronger dissipation. 
This means that within the hydrodynamic description there are time scales set by both the 
wavelength and the cooling rate and these can be quite different. 

A second question is the isolation of these hydrodynamic modes from the rest of the 
spectrum. At weak dissipation, it is expected that this is the case since the hydrodynamic 
eigenvalues are all small and the non-hydro dynamic spectrum is expected to be of the order 
of the collision frequency. At strong dissipation, there is a hydrodynamic mode of the 
order of Co/2, which can be of the order of the collision frequency. It remains an open 
question regarding the size of this eigenvalue relative to the non-hydrodynamic spectrum. It 
is possible that the isolation of the hydrodynamic spectrum places some restriction on the 
degree of dissipation. However, the analysis based on the kinetic model of Sec. EH suggests 
this may not be the case. For the kinetic model, all of the hydrodynamic spectra remain 
isolated from all of the non-hydrodynamic spectra for any degree of dissipation. In this case, 
the d + 2 hydrodynamic modes define a dominant set on a sufficiently large time scale, even 
when the separation of times among the hydrodynamic modes is significant. 

The idealized model for a granular gas discussed here can be made more realistic by con- 
sidering a more complex binary collision rule. In particular, a velocity dependent restitution 
coefficient (which approaches unity as the relative velocity goes to zero) and tangential fric- 
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tion are two additional qualitative features of real granular fluids. However, the general 
form of the hydrodynamic equations (macroscopic balance equations) is unchanged in that 
case; only the detailed values of the transport coefficients, pressure, and cooling rates differ 
from those of the present model. Similarly, the collision operator of the Boltzmann equation 
becomes more complex but its properties relevant for the macroscopic balance equations are 
the same. Most of the analysis given here depends only on those properties rather than the 
detailed form of the collision operator. Consequently, it is expected that the implications of 
our simple model extend to more realistic models as well. 

The analysis of the linear Boltzmann collision operator for normal, elastic gases is quite 
complete [lO|. It is hoped that the beginning provided here for granular gases will provoke 
the new mathematical analysis required for this well. 
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APPENDIX A: CHAPMAN-ENSKOG RESULTS 

The Chapman- Enskog method constructs a solution to the Boltzmann equation whose 
space and time dependence occurs entirely through the hydrodynamic fields and their gra- 
dients. For small spatial gradients, i.e. small relative variation of the hydrodynamic fields 
over a mean free path, the solution reads (we use the same notation as in ref. P|) 

f{r,v,t) = fl^J^{r,v,t)+A-V\nT + B-V\nn + C: Vu + MV^T + MW^n. (Al) 

Consider the contributions to the heat fiux and the momentum fiux from the terms of first 
order in the gradients of the above expression. Since these fiuxes appear in the macroscopic 
balance equations under a gradient, the resulting contributions are of second order in the 
gradients (Navier-Stokes order). For consistency, the cooling rate which does not occur 
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under a gradient in the balance equations must be calculated to second order. Therefore, 
the last two terms on the right hand side of Eq. ()A1|) give contributions to the cooling rate 
at this order, but lead to higher order (Burnett) terms in the heat and momentum fluxes. 
Additional non-linear terms of second order in the gradients coming from the cooling rate 
have been omitted in Eq. ()A1|) as only the linear hydrodynamic equations are considered 
here. The reference distribution function fj^J^{r,v,t) has the same functional form as the 
distribution function of the HCS discussed at the beginning of Sec. IIIIl but scaled with 
respect to the local exact hydrodynamic fields at time t for the generally non-uniform, non- 
stationary state. Therefore, it can be considered as the local HCS distribution function. 

The dissipative part of the pressure tensor and the heat flux are given by Eqs. and 
(jl2p . respectively. The transport coefficients in these equations are determined from the 
functions A, B, and C appearing in Eq. ()A1I ) through 



The functions A{v, n, T), B{v, n, T), and C{v, n, t) are solutions of the integral equations 





(A2) 



where the functions D{v) and S{v) are defined by 




(A3) 




(A4) 




(A5) 




(A6) 



with the definitions 



L'^'\t)X{v) = -J[fitx]-J[X,ffJ^], 



(A7) 
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d / ,(o)\ I d ( (0) 



Moreover, C,^^^^ is the local form of C/ics- 

The transport coefficients arising from the cooling rate are best described by first making 
the expression for the cooling rate more explicit. In general, it is given by a bilinear functional 
of the distribution function, 

C = C [/, /] , (A9) 

where 

= (1 - «') ^^^yr J d-"] dv,g^X{v)Y{v,) = C[Y,X]. (AlO) 



The linear contributions from the cooling rate to second order order in the gradients are 
then 

Cf = CiV2r + C2VV (All) 

with 

Ci = 2C[A^,/S], C2 = 2C[Ar,/(2]. (A12) 
The integral equations for Ai and A/" are found to be 

/ (0)^ 3 (0) (o)\ ^ _ df^ I 2k df^ J_ \ 

V ^''^^ dT 2^^^- + ^ )^~-^^' dT [dn dT ^ dT-^ ""J ' ^^^'^^ 

(0) rrAo) 



The above results can be easily transformed to be expressed in terms of the reduced units 
and quantities introduced in Sees. ITTland inTl Dimensionless transport coefficients are again 
defined by Eqs. fOHjl and (j^ . but replacing Uhcs and Thcs by their local values n and T. Of 
course, now it is £ = (na'''^^)^^ and Vhcsit) is replaced by v{t) = (2T/m)^/^. In this way, it 
is obtained: 

^* ^ ^2 _^ ^ _ 2 5Z (^^i' '^^i) ' ^ ^ 5Z '^»*) ' 
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d? 

with the dimensionless functions D* and S* given by 



D* = v*v* - ^ I, S* ={v*'^ - ] V*. (A16) 



Moreover, A*, B*, and C* are now defined through the equations 



C* -^]A* = A*, (A17) 



C*B* = B* + QA*, (A18) 



£* + ^ ) C = G*, (A19) 



where the operator C* is defined in Eq. (|^ and 

A* = v*^2(0) - *3(0), B* = 2v* -^s{0), 



G* = v**3(0) - ^v* ■ *3(0). (A20) 



The expressions of the reduced transport coefficients Q and Q are: 



where 



Cr = 3fa,A<*Y C2 = 3(«.^*)> (A21) 



2r(^) 

The integral equations obeyed hj M.* and A/"* are 



(1 - a2W(rf-i)/2 , 
= ^ L,, / dvl <p{vl)g*'. (A22) 



C* -k.]M* = (Cr - K*) ^2(0) + ^ ^* • v% (A23) 
2 / a 



£*Ar* = QM* + (C2* - /i*) ^2(0) + ^ i3* ■ t;*. (A24) 

For later use, it is convenient to elaborate more the above expression for (*. By con- 
struction, the velocity integrals of Ai times 1, v, and v"^ vanish. This is equivalent to say 
that is orthogonal to the set of functions $j defined in Eq. fj^ . and in particular it 
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verifies M.* = Q2M.*, with Q2 = I — V2, Vi being the projection operator defined in Eq. 
Then, acting with Q2 on both sides of Eq. ()A23|) it is obtained: 



M* = Q2M* = 

a 



Q2 ( /:* - f 



Q2A* ■ V* 



Q2 ( /:* - f 



For the same reason, Eq. (fMTjl yields 
A* = Q2A* = 

Q2 (^jC* 
Q2 (c* 



-1 



Q2A* 



2 

Co 



-1 



Q2 u* 



Q2V*^2(0) - 



Q2«*^2(0) + -*3(0). 

SO 



Co 



n -1 



*3(0) 



Substitution of this expression into Eq. ()A25j) after some algebra leads to 



M* = - 
a 



Q2 ( /:* - f 



Q2V* 



Q2 ( /:* - f 



and use of this into Eq. ()A21|) gives the result 



cr = ^ + i(<.,A^:). 



with 



Q2 ( - f 



Q2V* 



Q2 ( - f 



(A25) 



(A26) 



Q2t;*^2(0) - -2^1(0), (A27) 
Co 



(A28) 



Q2«*^2(0). (A29) 



Although the expression for ^2 can be written in a similar way, it will not be needed here. 



APPENDIX B: ADJOINT LINEAR OPERATOR AND BIORTHOGONAL SET 



The adjoint for C*, is defined as usual by 

(X, C*Y) = {C*^X, Y) 



(Bl) 



for arbitrary X{v*) and Y{v*) belonging to the Hilbert space. The explicit form of C* is 
given in Eq. (jSTj). From it, and using the above definition, the expression for C*'^ is easily 
found, 



C*^X{v*) = - J dvl4>{vl) j dae{a ■ g)a ■ g[X{v*") + X{v;") - X{v*) - X{vl)] 



2 dv* ^ 



(B2) 
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where v*" and t;*" are the postcoUisional velocities corresponding to v* and fjf, 

< = ^1* + ^ ■ 9*)^- (B3) 



Equation ()B2|) gives immediately 

r^l = 0, = -y^*, (B4) 

so that 1 and v* are eigenfunctions of the adjoint operator with eigenvalues and —Q/2, 
respectively. 

However, the kinetic energy is not an eigenfunction of the adjoint operator. Direct cal- 
culation gives 

where a{v) is given by Eq. ()A22j) . Nevertheless, a biorthogonal set can be constructed from 
1, V*, and f*^ and it is given in Eq. (j3Hj) . As noted in the main text, the choice of this set 
is not unique. The conditions of biorthogonality on $2 are 

<!>2,^m) = (v*-^,l)=l, (B6) 



$2,^1(0) =($2,1) -1 = 0, (B7) 



"^2,*3(0))= =0. (B8) 



dv* 

A sufficient condition to guarantee that the above relations are verified, is that $2 have the 
form 

<^2{v*)=A + Bb{v*), (B9) 

where b{v*) is an arbitrary scalar function of f *, and the constants A and B are determined 
from ^ 

^ = l-(l,6)i?, B=(^v*.^,l^ . (BIO) 

Substitution of these expressions into Eq. ()B8|) yields 



^,iv*) = i + [b{vn-ii,b)]{v*--^,i] . (Bii) 
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The optimal choice for b{v*) would be that implying that ^2{v*) is an eigenfunction of C*^ 
corresponding to the eigenvalue Co/2- This is accomplished if b{v*) is the solution to 

£%{v*) = ^ [b{v*) + B-^] . (B12) 

The solution to this equation, if it exists, has not yet been found. 



APPENDIX C: EVALUATION OF THE PERTURBATION THEORY RESULTS 

The expansion of the hydrodynamic eigenvalues of the linearized Boltzmann equation for 
small k is given in Sec. |V|with the result 

A,(k) = A,(0) + A;2Af + ..., (CI) 



where 



and 



{A.(0)} = ^0,|,-|^, (C2) 



k ■ V* {Q, [C* - K Q^k ■ v*^, (0)^ + ($„ r Q,^f ) . (C3) 



The eigenvalue —Q/2 is d-io\d degenerated and the convenient choice for the lowest order 
eigenfunctions has been discussed in Sec. IIVI The formal expression for the second order 
eigenfunctions Qi'^l is given in Eq. ()70|). For i 7^ 2, the second term on the right hand side 
of Eq. ()C3|) vanishes since $i is an eigenvector in those cases and, therefore. 



(<l>i, £*Qi^f ^) = (^C*^<l>i, Q.^f ) oc ($„ Qi^P) = 0. (C4) 
Then, Eq. ()C3|1 can be rewritten as 

where it has been used that 

Pifc ■ v*^i{S)) = ^i(O) k ■ i;*^i(0)^ = 0. (C6) 

This follows since $j and \l'i(0) have the same parity with respect to reflections of v* . 
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For the first eigenvalue it is 



(C7) 



1 

So V /So 
In the first transformation, the property ViC*X = 0, for arbitrary X{v*), has been em- 
ployed. Next, consider the eigenvalue associated to the longitudinal component of ^'3(0) 
that we will denote by An , 



.(2) 



kv' 



Co 



n -1 



k ■ v*k 



din 4 
dv* 



*2 



Q3 + 



So 



1 -1 



Q3 ( + f 



1 -1 



d+ 1 

^1(0) + ^- ^2(0) 



*2 

7; 



d Co 



(C8) 



where it has been taken into account that {C* + Q/2)=0. Using Eq. ()A19|) in the first 
term on the right hand side of the above equation and evaluating explicitly the other two 
ones one gets 

(C9) 



The last equality follows from the fact that the expression of rj* given in Eq. ()A15|) is 
equivalent to 



d-2 



6, 



(CIO) 



because of the symmetry of the tensors D and C. The calculation of the shear modes 
(2) 

eigenvalues is straightforward: 



A 



(2) 



So 



3,±i 



2 / 



(Cll) 
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(2) 

Finally, the evaluation of is somewhat more complicated. The contributions from each 



of the two terms in Eq. ()C5|) will be computed separately. The first one is given by 



Q2 - f 



-1 



fc-«*vi/2(o)) = (<<i>2,(/:-|) v:^,io)], (C12) 



that is equivalent to 



= I (SI A*) - ^ {s:, ^,M) - ^ ^ ( <' ^^2(0) 



(C13) 



The analysis of the second term on the right hand side of Eq. ()C5|) is carried out in an 
analogous way, 

\ {C*v*\ Q.M/f) = -| (^*^ Q^v&r) + i (a, . (CM) 



It is 



^ (^*^ Q^^?^) = I ($1, 22^1^?^ 

1 _ _ , 1 



-($3„<v[;2(0)) = --. (C15) 

S,0 So 



In the second equality above, the explicit expression of \E'2 given in Eq. (fTOj) as well as the 
properties of C* have been used. The second term on the right hand side of Eq. ()C14|) can 
be rewritten as 

1 / - .,(2)^ 1 



where Ml is defined in Eq. flX29)l . Substitution of Eqs. (IHTsl) and KJIM into Eq. dnTijl 
yields 

^ (/:v^ Q^vi^f ) = - ^ = -cr, (ci?) 

because of Eq. (1X281) . Then, use of Eqs. (inT3)l and (inT7|l into Eq. JCSl) gives the final 
expression for Aa^"*, 

Af = - C - ^7^. (C18) 

Comparison of the results for the second order eigenvalues obtained in this Appendix with 
those reported in Sec. ITTl shows that both agree, with the transport coefficients given by the 
same expressions in both cases. 
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APPENDIX D: LINEARIZATION OF THE MODEL KINETIC EQUATION 



Solutions to the model kinetic equation (f?^ are sought of the form (jH^ and, consistently, 
9 = Qhcs + 5g and v = v^cs + ^i^i retaining only terms linear in A, 



d 



Use has been made of the HCS model equation 

dfhcs Chcs ^ ( f \ ( f \ 

^ ' y^Jhcsj ^hcs yjhcs Qhcsj ■ 

Next, from the choice of i> as corresponding to hard sphere behavior. 



/ 5n 26T 
oy = Vhcs I \- TfT- 

\^hcs ^hcs 

Similarly, the definition of g in Eq. ()73p gives 



r- QhcS r- dghcs ;m 9ghcs 

dg = 6n + — — 6T - — — ■ u. 

n-hcs oThcs ov 

The linear model kinetic equation ()D1|) then becomes 



^ + v-V)(A,,A) 



6n 

^hcs 



6T 



2T, 



has 



^hcs i^fhcs 9hcs) 



I £ A 9hcs r- dghcs ;m , ^dhcs 

-J^hcs Acs A dn - — — dT + — — 

rihcs oThcs ov 



u 



It is convenient to eliminate ghcs in this result using the HCS equation ()D2|) to get 

— + V ■ V (Acs A) = -VhcsJhcs^ + I Z ^ T;?^ ) — — ■ V"Jhcs) 



dt 



rthc 



2T, 



hcs 



2 dv 



6n 
+ 

TT'hcs 

+6T 



hcs fhcs + ^7;~7^ ■ iyfhcs) 



2 dv 

dfhcs Chcs d 
^hcs-^ r 



'dTj 



hcs 



V- 



df 



hcs 



— U ■ 



d_ 
dv 



^hcsfhcs 



2 dv y dThcs 

Chcs d 



2 dv 



(Vfhcs) 



(Dl) 



(D2) 



(D3) 



(D4) 



(D5) 



(D6) 



The right hand side of this expression can be put in a more transparent form by noting that 



6n 6T Sn u\\ Ui 



T^hcs hcs TT'hcs '^hcs '^hc 



(D7) 
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where the scalar product {X,Y) is defined in Eq. (pUjl and are the biorthogonal set 
given in Eq. (j3H|) . Moreover, the relation 



T; 



d\nfi 



hcs 



1 d 



hcs 



[V6 



(D8) 



dThcs 2(1) dv 

allows to write the dependence on fh^s on the right hand side in terms of the lowest order 
hydrodynamic eigenfunctions {\l'i(0)} of the linearized Boltzmann operator given in Eq. 
()54|) . Then, Eq. ()D6|) can be rewritten as 

d 



dt 



hcs 



(hcs d 

Vhcs{t) 



hcs ' 



A-5^vl/^(0) ($„A) 

i 

fhcsVj^'^M {'^^,^) 

i 



(D9) 



Here {Aj(0)} are the eigenvalues of Eq. 

Finally, introducing the dimensionless variables of Sees. ITTland inTl and the single Fourier 
component of Eq. ()35|) gives the linear kinetic equation for the model 



with 



{ds + tk-y^*+C*JA{k,v*,s) = 0, 



2 " dv* 

The projection operators Vi are defined in Eq. (jU^ . and 



It is directly verified that 



i 

-1 d 



dv* 



v*4>Q = 0, 



(DIO) 



(Dll) 



(D12) 



(D13) 



so that can be written 



V 



Q. 



2 ^ dv 

This is the expression for the generator of the linear dynamics used in main the text. 



(D14) 
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